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I. INTRODUCTION 

Since the pioneering work of the Anderson 1 we know 
that disorder localizes electrons. At zero temperature 
any small amount of the disorder localizes all electrons in 
the one-dimensional systems 2 . In higher dimension (d > 
2 for systems with orthogonal symmetry) weak disorder 
does not destroy the metallic regime. Only when the 
strength of the disorder increases over a certain specific 
value (critical disorder), the electron becomes localized. 
This phenomenon - transition from the metallic to the 
insulating regime due to an increase of disorder - is called 
Anderson transition 3 . 

Scaling theory of the Anderson transition uses the con- 
ductance g 4 as the order parameter. It is supposed 5-7 
that the system size dependence of the conductance is 
determined only by the value of the conductance itself: 



d In g 
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(1) 



where 13(g) is an analytical function of g. (3 is posi- 
tive (negative) in the metallic (localized) regimes, respec- 
tively. For dimension d > 2 the function (3(g) changes its 
sign, being positive for j > 1 and negative in the limit 
g — ► — co. There is an unstable fixed point g c defined as 
the solution of (3(g = g c ) = 0. System-size independent 
critical conductance g c represents the critical point of the 
Anderson transition. 

Relation (1) contains no information about the micro- 
scopic structure of the model. This means that the An- 
derson transition is universal. The form of the /3 function 
is determined only by the physical symmetry and dimen- 
sion of the system 8 . 

Soon after the formulation of the scaling theory of lo- 
calization it became clear that the conductance g is not a 
self-averaged quantity. Reproducible fluctuations of the 
conductance were found both in the metallic and in the 
insulating regimes 9 ' 10 . The knowledge of the mean value 
(g) is therefore not sufficient for complete description of 
the transport properties. One has to deal with the con- 
ductance distribution P(g) 11,12 or, equivalently, with all 
cummulants of the conductance. This is easier in the 
metallic regime, where P(g) is Gaussian and the conduc- 
tance fluctuations are universal 13-21 and independent on 
the value of the mean conductance and/or the system 
size. The width of the distribution depends only on the 
dimension, physical symmetry of the system 13 and on the 
boundary conditions 22,23 . In the insulator, conductance 
wildly fluctuates within the ensemble of macroscopically 
equivalent ensembles. It is the logarithm of the conduc- 



tance which is distributed normally in the the limit of 
large system size 11 ' 12 ' 17,24 . 

At present we have no complete analytical theory 
able to describe the conductance statistics at the criti- 
cal point. Analytical results are known only for systems 
of dimension d = 2+e (e <C l) 12 > 25 > 26 . Expression for con- 
ductance cummulants 25 enabled to estimate the shape of 
the critical conductance distribution P c (g) 26 . P c (g) is 
system size independent and Gaussian in the neighbor 
of the mean value (g) ~ (T 1 . The distribution possesses 
power law tail P c (g) ~ <? -1-2 ^ e for g — > oo, and the 
delta-function peak 8(g). These results can not be ap- 
plied to three dimensional (3D) system (e = 1) where 
(g) ~ l 27 . Numerical simulations are therefore crucial 
in 3D system. The first systematic numerical analysis of 
the conductance statistics in 3D was done in 27 and was 
followed by a scries of papers 28-34 . In 2D systems, crit- 
ical conductance distribution was numerically studied in 
the regime of quantum Hall effect 35 and in systems with 
the spin-orbit interaction 28,36,37 . 

This paper reviews our recent numerical data for the 
conductance distribution. We address the question of the 
shape, universality and the scaling of the critical conduc- 
tance distribution. 

In numerical calculation of the conductance we sup- 
pose that two opposite sites of the sample are connected 
to semi-infinite perfect leads and use the multichannel 
Landauer formula 38 which relates the conductance g (in 
units 2e 2 /h) to the transmission matrix t: 
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In (2) we introduced the variables Zj, i — 1,2... N open , 
(zi < Z2 < ■ ■ •) which parametrize the eigenvalues of the 
matrix tH. In the limit L z >> L, Zi converges to 2L Z /Xi 
where is the ith localization length of the quasi-one 
dimensional (quasi-lD) system 17,27 
of open channels. 

Owing to relation (2), the analysis of the conductance 
can be reduced to the calculation of the eigenvalues of 
the matrix tH. The general formula for t was derived 

m 39,40 

According to the scaling theory, critical exponents of 
the Anderson transition as well as conductance statistics 
do not depend on the microscopic details of the model. In 
numerical simulations, we consider the Anderson Hamil- 
tonian 
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FIG. 1. Left The L - dependence of the mean conductance 
(g) for different values of the disorder in the critical region. 
Data show that 16.45 < W c < 16.55. More exact estima- 
tion of the critical disorder was done in 47 . Right The un- 
ambiguous dependence var g vs (g) in the critical region of 
the metal-insulator transition. This agrees with the single 
parameter scaling theory. 

In (3) n counts the sites of the d-dimensional lattice and 
[nn'\ means nearest neighbor sites. The hopping term r 
equals to unity for orthogonal systems, unless the ani- 
zotropy is considered. When spin orbit scattering is con- 
sidered, t is a 2 x 2 matrix 41 ' 42 . Parameter W measures 
the strength of the disorder. For box distribution of ran- 
dom energies e„ (P(e) = 1 for |e| < 1/2 and P(e) = 
otherwise), the 3D Anderson model (3) exhibits Ander- 
son transition at W r = 16.5. 



numerically for the 3D Anderson model 47 . Statistical en- 
sembles of more than -/V sta t > 10 6 cubes of the size from 
4 3 to 18 3 were collected for various values of the disorder 
W. Fig. 1 shows typical data for the system size de- 
pendence of the mean conductance. In agreement with 
(1), (g) increases (decreases) with the system size in the 
metallic (localized) regime. By the use of the general 
fitting procedure 46 , the critical disorder and the critical 
exponent v 1.57 were obtained. 

Data in fig. 1 confirm that the variance, var g — 
(g ) — (g) 2 is an unambiguous function of the mean (g) 
in the critical regime. This supports, but still does not 
prove the single parameter scaling theory. General proof 
of the theory requires verification of the single parame- 
ter scaling of all conductance cummulants. This is nu- 
merically impossible since higher cummulants are fully 
determined by rare events with very large values of the 
conductance. 



B. Scaling of the conductance distribution 

As higher cummulants are not treatable numerically, 
we test the scaling of the conductance distribution by 
the analysis of the scaling of percentiles g a 34 - Percentile 
g a is defined as 



II. FINITE SIZE SCALING 

As it was discussed in the Section I, the conductance 
g is not the self-averaged quantity. To avoid statisti- 
cal fluctuations, quasi-lD systems were introduced 43 and 
the localization length 7 is calculated instead of the con- 
ductance. In the limit L z /L — > 00 (L z and L are the 
length and the width of the system, respectively) 7 is a 
self-averaged quantity which converges to its mean value. 
Finite size scaling is then used 44 for the analysis of the 
disorder and the system width dependence of 7. It is 
assumed that the variable A = 7/L is a function of only 
one parameter: A(L,W) = A(L/£(W)). Here £ = £(W) 
is the correlation length which diverges in the vicinity of 
the critical point £(W) = \W— W c \~ u . Critical exponents 
v and s = (d— 2)i> characterize the critical behavior of the 
localization length and of the conductance, respectively 5 . 

Finite size scaling analysis of the quasi-lD data en- 
abled to test the universality of the Anderson transition 45 
and provided us with the more accurate estimation of the 
critical exponent v sa 1.57 46 . 



A. Scaling of the mean conductance 

Verification of the scaling theory of the localization 
requires the proof of the universal scaling of the mean 
conductance and of the entire conductance distribution 
in the critical regime. Single parameter scaling of vari- 
ous mean values, (g), exp(lng), and l/(g _1 ), was proved 



P(g)dg. 



(4) 



Owing to (4), the probability to find g < g a equals to a. 
Of course, the percentile g a is a function of disorder and 
system size: g a — g a (L, W). Single parameter scaling of 
percentiles has been proved for several values of a 34 . 
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FIG. 2. Left Critical conductance distribution P c {g) of the 
3D Anderson model, obtained for statistical ensembles of 
A^stat = 10 6 samples of the size 10 3 - 18 3 . P c (g) is sys- 
tem size independent and decreases faster than exponentially 
when g > 1. The main properties of the critical distribu- 
tion are shown in insets: Left inset shows in details the small 
- g behavior and proves that the distribution decreases to 
zero as g — > + . Right inset shows the discontinuity of the 
derivation dP(g)/dg at g = 1. Right figure presents the dis- 
tribution P c (ln(?) at the critical point for the three dimen- 
sional (W c = 16.5, open symbols) and four dimensional (4D) 
(Wc = 34.3, full symbols) systems of various system size. 
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FIG. 3. Conductance distribution for ensembles of cubes 6 3 
(dashed lines) and 12 3 (solid line): (a) The metallic regime 
W — 15.4 (the critical disorder is W c = 16.5). Distribu- 
tion P(g) moves toward higher values of the conductance as 
the system size increases and P(g) becomes Gaussian when 
L — > oo. (b) Insulator W = 17.6 Mean conductance decreases 
as the system size increases and the distribution of the lng 
becomes Gaussian in the limit L —* oo. 
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FIG. 4. Comparison of the critical conductance distribu- 
tion of two different 2D models with spin-orbit coupling: 
full symbols Ando model 42 open symbols Evangelou-Ziman 
model 41 . Squares 40 x 40 and 80 x 80 were simulated to prove 
the system size independence of the critical distribution. Note 
also the non- analytic behavior for g = 1. 



Suppose that g a and gp (a < (3) obey the single pa- 
rameter scaling. Then g 1 (a < 7 < (3) scales, too. There- 
fore, in contrast to the analysis of the conductance cum- 
mulants, it is enough to analyze only a few percentiles. 
Next, if g a and gp scale, then the difference gp — g a scales. 
Scaling of percentiles assures thus the scaling of the en- 
tire conductance distribution. Of course, this analysis is 
not applicable to the limit a — > 1, because available sta- 
tistical ensembles are never big enough to provide us with 
sufficient information about the tail of the distribution. 



III. CRITICAL CONDUCTANCE DISTRIBUTION 

At present, we have no analytical description of the 
critical conductance distribution in 3D systems. Analyt- 
ical results were obtained only for the conductance cum- 
mulants in the dimension d = 2 + e close to the lower crit- 
ical dimension (e <C l) 25 . In spite of the non-universality 
of higher order conductance cummulants 

0' > = {~£~"- „> e -l (5) 

the critical distribution P c {g) was shown to be universal 
and L-independent in the limit L — > oo 26 . However, the- 
oretical analysis of the form of the critical distribution, 
is applicable only in the limit of very small e 27 . 

A. The form of the critical conductance distribution 

All what we know about the P c (g) in 3D is based 
on the numerical data. In Fig. 2 we present P c (g) for 
the 3D Anderson model. Data confirm that the critical 
conductance distribution is system size independent, as 
required 12 . The shape of P c {g) differs considerably from 
the conductance distribution in the metallic and in the 
insulating regimes. To explain the typical properties of 
the critical conductance distribution, we use our knowl- 
edge about statistical properties of parameters z (2). 



Decrease of P c (g) to zero when g — > is visible when 
very large statistical ensembles are studied (left inset of 
Fig. 2). Due to Eq. (2), small conductance means that Z\ 
is large. From numerical data we know that the distribu- 
tion P(z\) is similar to Wigner surmises and decreases as 
exp — z? for large z\. Consequently, lnP c (ln<7) decreases 
as — In g (right figure in Fig. 2), and guarantees that 
lim s ^ P c ( 5 )= 27 ' 30 . 

Large g behavior of P c {g) is determined by the chance 
that many parameters Zi are small. Statistical analysis 
of parameters zt showed that for i > 1 the distribution 
P(zi) is Gaussian with mean value (zi) cx yfi and var 
Z{ oc (z;)~ 130 . The probability to find a sample with 
small value of higher Zi is therefore very small: P c (g ~ i) 
decreases as exp— i 3 / 230 . Fig. 2 indeed shows very fast 
decrease of the probability P c (g) for g > 1. The chance 
to have g > 1 is only 3%. Probability to find large values 
of g drastically decreases: we found that in the ensemble 
of 10 7 samples (L = 10) only 470 samples have g > 2 and 
only one sample has g > 3 33 . The analysis of the con- 
tribution of the first two channels is therefore sufficient 
for the understanding of the qualitative properties of the 
critical distribution. 

Numerical data also show that the critical distribution 
is non-analytical at g = 1. Right inset of Fig. 2 shows the 
discontinuity of the first derivative dP c (g)/dg. The same 
non-analyticity was found in the 4D systems (right figure 
2), in the unitary 29 and symplectic (Fig. 5) systems, and 
also in the weakly disordered quasi-lD systems 48 ' 49 . 

Present description of the critical distribution is based 
on the analysis of statistical properties of parameters z. 
It is applicable to any system, for which the mean values 
of parameters z are of the order of unity. Then only a few 
(two, or three) channels contribute to the conductance. 
This analysis is, however, not applicable to systems close 
to the lower critical dimension d = 2 + e. Here, the 
mean conductance (g) ~ e~ x , which means that number 
of channels which contribute to the conductance, is large, 
~ e _1 . It is therefore no surprise that the critical con- 
ductance distribution found in 26 differs from that shown 
in fig. 2. 
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FIG. 5. Critical conductance distribution for bifractal lat- 
tices. Legend gives spectral dimension. 5th generation of 
fractals was used to calculate the distribution. Bifractal is 
linear along the z axis and fractal in a plane perpendicular 
to z. Right figure shows the structure of the fractal (the 3rd 
generation) cross section. All fractals have the same fractal 
dimension df — In 3/ In 2. Spectral dimension d s — In 9/ In 5 
for fractals A and B, and = In 9/ In 6 for fractal C. Note that 
although bifractals A and B have the same both fractal and 
spectral dimensions, they posses different critical conductance 
distribution. P c {g) depends on the lattice topology 51 

B. Conductance distribution in the critical regime 

Fig. 3 shows that the typical properties of the crit- 
ical conductance distribution hold also for P(g) in the 
neighbor of the critical point. This is because the statis- 
tical properties of parameters z depend continuously on 
the disorder in the critical regime. Critical properties of 
P(g) survive until the system size L exceeds the correla- 
tion length £. Only when («L the distribution typical 
for the metallic or the insulating regime can be observed. 

C. Universality 

Single parameter scaling theory of localization sup- 
poses that the critical conductance distribution is uni- 
versal. Its form does not depend on the microscopic de- 
tails of the model. As an example, we present in Fig. 
4 the critical conductance distribution for two 2D mod- 
els with spin-orbit scattering: Evangelou-Ziman model 
and Ando model 42 . In spite of the different microscopic 
definition of both models, P c (g) is universal 28 . Univer- 
sality of P c {g) with respect to various distributions of 
random energies e„ was confirmed in 28 and 32 . 

As was shown already in fig. 2, the shape of P c (g) de- 
pends on the dimension of the model. Ref. 29 confirmed 
that also the physical symmetry influences the form of 
P c (g)- Less expectable was the observation 22 ' 31 that 
P c (g) as well as the spectral statistics 50 depend on the 
boundary conditions in the transversal direction. Never- 
theless, this is consistent with the original definition of 
the conductance as a measure of the sensitivity of the en- 
ergy spectrum of the system to the change of the bound- 
ary conditions 4 . 

For completeness, let us note that the critical con- 
ductance distribution depends also on other parameters 
of the model: lattice topology 51 , anizotropy 37 and, of 



course, on the length of the system. We believe that these 
non-universalities could be compensated by the change of 
another model parameter (see, for instance 37 , where the 
anizotropy is compensated by the length of the system). 

D. Dimension dependence 

Right figure in Fig. 2 compares the critical con- 
ductance distribution for 3D and 4D cubes. As sup- 
posed, the maximum of P c {g) for 4D is shifted toward 
smaller conductances, because the critical disorder in- 
creases as the spatial dimension increases and higher dis- 
order means lower mean conductance 5,26 . Qualitatively, 
however, both distributions are very similar: P c (ln<7) de- 
creases as exp[— In 2 g] for In g — > — oo and possesses the 
non-analyticity at ln g = 0. This similarity is not surpris- 
ing, because the form of the distribution is determined 
mostly by the statistics of Z\ and Z2, which are qualita- 
tively similar in 3D and 4D 30 . 

Surprisingly, the relation (5) seems to hold also for 
e = 1 and 2, at least for the first two cummulants. We 
obtained numerically that (g) = 0.285 for 3D and 0.135 
for 4D, so that (g)$D ~ 2(3)40. For the second cummu- 
lants we found var g^D w varg4D ~ 0.17 . 

More interesting is the investigation of the P c {g) in 
the systems of dimension 2 + e 51 . As we are not able to 
create the d-dimensional hyper-cubes with non-integer d 
in computers, we simulated the transport on bifractal 
latices 52,51 . Bifractals are linear along the propagation 
direction and possess the fractal lattice in the cross sec- 
tion (fig. 5). We proved that the critical exponent v de- 
pends only on the spectral dimension of the lattice. Mean 
conductance, var g and the critical distribution P c {g) de- 
pend, however, on the lattice topology. For instance, fig. 
5 shows that bifractals A and B have the same spectral 
dimension, but different critical distribution. This is the 
reason why obtained data can not be used for the verifi- 
cation of relations (5). 

In Fig. 5 we present P c (g) for three different bifractals. 
As expected, (g) increases and the distribution converges 
to Gaussian when e — > 0. However, we found neither the 
S - function peak at g = nor the power-law tail of the 
distribution for g ^> (g) , predicted by the theory 26 . 

IV. CONDUCTANCE DISTRIBUTION IN 
NON-CRITICAL REGIME 

Although we have no analytical theory of the con- 
ductance statistics in the critical regime, we can learn 
some typical properties of the conductance distribution 
from the analysis of the quasi-lD weakly disordered 
systems 48 . Starting from the Dorokhov - Mello-Pereyra- 
Kumar equation 15 for the probability distribution of pa- 
rameters 2, the conductance distribution P(g) can be 
calculated. P(g) depends on the length of the system. 
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FIG. 6. Comparison of the conductance distribution of the 
3D (dashed lines) and quasi- ID (solid line), (a) Metallic 
((g) > 1) (b) critical ((g) « 1) and (c) localized ((g) < 1) 
regimes are shown. In the quasi-lD systems, the strength of 
the disorder is fixed to W — 4 and the length of the system 
is tuned to obtain the same mean conductance as in the 3D 
system. 
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FIG. 7. Distribution of — In g for strongly localized regime. 
Note that the distribution is not symmetric. It possesses the 
long tail for small values of the argument, and decreases much 
faster on the opposite side. Solid lines are a Gaussian distribu- 
tion with the mean and variance found from numerical data. 
Size of the cubes is L = 10. Left W = 38. Note that the 
sharp decrease of the distribution at \ng — 0, discussed in 
Sect. 3.1, is still present 48 . Right W = 60. 

For short samples, P(g) is Gaussian 17 (see also fig. 6a). 
When the length of the system increases, the conduc- 
tance decreases and the system passes from the metallic 
regime into the localized one. For the intermediate sys- 
tem length, where (g) w 1, P(g) is expected to be quali- 
tatively similar to the critical conductance distribution 48 . 
This is shown in Fig 6, where we compare the conduc- 
tance distribution of quasi-lD and 3D systems of the 
same mean conductance. 

Fig. 6 also shows that the analogy between 3D and 
quasi- ID systems is not exact. To understand quantita- 
tive differences between the quasi- ID and 3D systems, we 
analyzed the spectrum of parameters z. In the metallic 
regime, the difference between 3D and quasi-lD systems 
is only in the value of var g 13 . P(g) is Gaussian and the 



dently on the dimension of the system 54 . In contrast to 
the metallic regime, the spectrum of z becomes dimen- 
sion dependent in the critical region. Both quasi-lD 53 ' 54 
and 3D-dimensional 49 numerical studies confirmed that 
at the critical point 



cx i 



(6) 



(d > 2). Owing to (6), the difference A = (z$ — z\) 
is smaller that (z%) in 3D, while it equals to (z%) in the 
quasi- ID system. The contribution of the second channel 
is therefore more important in 3D than in the quasi-lD. 
This explains longer tail of the distribution in the 3D 
system (fig. 6b). 

It is commonly believed that the distribution P(\ng) 
is Gaussian in the insulating phase, independently on 
the dimension of the system. This is, however, not true. 
The spectrum of z depends namely on the dimension of 
the system also in the localized regime. For the 3D sys- 
tems it was proved numerically 49,54 that the difference 
A = (z% — Zi) is constant, independent on the disorder 
and on the system size. Therefore the second channel 
influences always the statistics of the conductance. More 
than its contribution to the value of the conductance it 
is important that constant value of A prevents the distri- 
bution P(z\) to develop into the Gaussian form. While 
the values z\ <C (zi) are still possible, the probability 
to find systems with much higher values Z\ 3> (z\) is 
strongly suppressed. The distribution P(z\) is therefore 
not symmetric. The same is true for the distribution 
P(—\ng) (Fig. 7) which possesses a long tail for small 
values of \ \xig\ and decreases much faster than Gaussian 
for \ hxg\ — > oo. 

Note that in weakly disordered quasi- ID systems A ~ 
(zi). The distance (z2 — z\) is much larger than the width 
of the distribution P(z\). Higher channels therefore do 
not influence the distribution of z\ and P In g is Gaussian. 



V. CONCLUSION 

We reviewed recent progress in numerical studies of 
the statistics of the conductance in the critical regime. 
Numerical analysis confirms that the conductance dis- 
tribution in 3D Anderson model obeys single parameter 
scaling. Analysis of the statistics of the eigenvalues of the 
transmission matrix enables us to understand the main 
features of the conductance distribution in the critical 
regime. Although we still have no analytical description 
of the conductance statistics in the critical regime, we 
hope that results of numerical experiments will inspire 
theoreticians to formulate the general analytical theory 
of the Anderson transition. 



spectrum of parameters z is linear 



indepen- 



5 



1 P. W. Anderson: Phys. Rev. 109, 1492 (1958) 

2 N. F. Mott, W. D. Twose: Adv. Phys. 10, 107 (1961) 

3 N. F. Mott, E. A. Davis: Electron processes in non- 
crystalline materials, Clarendon Press, Oxford 1979 

4 J. T. Edwards, D. J. Thouless: J. Phys. C 5, 807 (1972); 
D. C. Licciardello, D. J. Thouless: J. Phys. C 8, 4157 (1975) 

5 F. Wcgner: Z. Phys. B 25, 327 (1976) 

6 E. Abrahams, P. W. Anderson, D. C. Licciardello, 
T. V. Ramakrishnan: Phys. Rev. Lett. 42, 673 (1979) 

7 L. P. Gorkov, A. I. Larkin, D. E. Khmelnitskii: JETP Lett. 
30, 228 (1979) 

8 B. Kramer, A. MacKinnon: Rep. Prog. Phys 
56, 1469 (1993) 

9 S. Washburn, R. A. Webb: Adv. Phys. 35, 375 (1986) 

10 A. B. Fowler, J. J. Wainer, R. A. Webb: IBM J. Res. Dev. 
32, 372 (1988) 

11 P. W. Anderson, D. J. Thouless, E. Abrahams, D. S. Fisher: 
Phys. Rev. B 22, 3519 (1980) 

12 B. Shapiro: Philos. Mag. B 56, 1032 (1987); Phys. Rev. 
Lett. 65 1510 (1990) 

13 P. A. Lee, A. D. Stone: Phys. Rev. Lett. 55, 1622 (1985); 
P. A. Lee, A. D. Stone, H. Fukuyama: Phys. Rev. B 
35, 1039 (1987) 

14 I. Imry: Europhys. Lett. 1, 249 (1986) 

15 O. N. Dorokhov: JETP Lett. 36, 318 (1982); P. A. Mello, 
P. Pereyra, N. Kumar: Ann. Phys. (NY) 181, 290 (1988) 

16 C. W. J. Beenakker: Rev. Mod. Phys. 69, 731 (1997) 

17 J.-L. Pichard: 'Random transfer matrix theory and con- 
ductance fluctuations'. In: Quantum Coherence in Meso- 
scopic Systems NATO ASI Ser. B: Physics vol. 254 ed. by 
B. Kramer (Plenum Press, 1991) 

18 P. M. Bell, A. MacKinnon: J. Phys.: Condens. Matt. 
6, 5423 (1994) 

19 N. Giordano: Phys. Rev. B 38, 4746 (1988) 

20 J.-L. Pichard, N. Zannon, I. Imry, A. D. Stone: J. Phys. 
France 51, 587 (1990) 

21 Y. Avishai, J.-L. Pichard, K. A. Muttalib: J. Phys. France 
3, 1387 (1991) 

22 D. Braun, E. Hofstetter, G. Montambaux, A. MacKinnon: 
Phys. Rev. B 64, 155107 (2001) 

23 M. Riihlander, P. Markos, C. M. Soukoulis: Phys. Rev. B 
64, 172202 (2001) 

24 P. Markos, B. Kramer: Ann. Physik (Leipzig) 2, 339 (1993) 

25 B. L. Altshuler, V. E. Kravtsov, I. V. Lerner: Sov. Phys. 
JETP 64, 1352 (1986); JETP Lett. 43, 441 (1986) 

26 A. Cohen, B. Shapiro: Int. J. Mod. Phys. B 6, 1243 (1992) 

27 P. Markos, B. Kramer: Philos. Mag. B 68, 357 (1993) 

28 P. Markos: Europhys. Lett. 26, 431 (1994) 

29 K. Slevin, T. Ohtsuki: Phys. Rev. Lett. 78, 4083 (1997); 

30 P. Markos: Phys. Rev. Lett. 83, 588 (1999) 

31 K. Slevin, T. Ohtsuki, T. Kawarabayashi: Phys. Rev. Lett. 
84, 3915 (2000) 

32 K. Slevin, T. Ohtsuki: Phys. Rev. B 63, 45108 (2001) 

33 M. Riihlander, P. Markos, C. M. Soukoulis: Phys. Rev. B 
64, 212202 (2001) 

34 K. Slevin, P. Markos, T. Ohtsuki: cond-mat/0208157 

35 Z. Wang, B. Jovanovic, D. H. Lee: Phys. Rev. Lett. 
77, 4426 (1996); V. Plerou, Z. Wang: Phys. Rev. B 
58, 1967 (1998); X. Wang, Q. Li, C. M. Soukoulis: Phys. 
Rev. B 58, R3576 (1998) 



36 P. Markos: J. Phys. I France 4, 551 (1994) 

37 M. Riihlander, P. Markos, C. M. Soukoulis: Phys. Rev. B 
64, 193103 (2001) 

38 R. Landauer: IBM J. Res. Dev. 1, 223 (1957); R. Landauer: 
Phil. Mag. 21, 683 (1970); E. N. Economou, C. M. Souk- 
oulis: Phys. Rev. Lett. 46, 618 (1981) 

39 T. Ando: Phys. Rev. B 44, 8017 (1991) 

40 J. B. Pendry, A. MacKinnon, P. J. Roberts: Proc. P. Soc. 
London A 437, 67 (1992) 

41 S. N. Evangelou, T. Ziman: J. Phys. C 20, L235 (1987) 

42 T. Ando: Phys. Rev. B 40, 5325 (1989) 

43 J.-L. Pichard, G. Sarma: J. Phys. C 14, L127 (1981); ibid, 
L167 

44 A. MacKinnon, B. Kramer: Phys. Rev. Lett. 
47, 1546 (1981) 

45 B. Bulka, B. Kramer, A. MacKinnon: Z. Phys. B: Condens. 
Matt. 60, 13 (1985); B. Kramer, K. Broderix, A. MacKin- 
non, M. Schreiber: Physica A 167, 163 (1990); A. MacK- 
innon: J. Phys.: Condens. Matt 6, 2511 (1994) 

46 K. Slevin, T. Ohtsuki: Phys. Rev. Lett. 82, 382 (1999) 

47 K. Slevin, P. Markos, T. Ohtsuki: Phys. Rev. Lett. 
86, 3594 (2001) 

48 K. A. Muttalib, P. Wolfle: Phys. Rev. Lett. 83, 3013 (1999); 
V. A. Gopar, K. A. Muttalib, P. Wolfle: cond- 
mat/0205643; K. A. Muttalib, P. Wolfle, A. Garci'a-Martm, 
V. A. Gopar: cond-mat/0210293 

49 P. Markos: Phys. Rev. B 65, 104207 (2002) 

50 D. Braun, G. Montambaux, M. Pascaud: Phys. Rev. Let. 
81, 1062 (1998) 

51 I. Travenec, P. Markos: Phys. Rev. B 65, 113109 (2002) 

52 M. Schreiber, H. Grufibach: Phys. Rev. Lett. 
76, 1687 (1996) 

53 P. Markos, M. Henneke: J. Phys.: Condens. Matt. 
6, L765 (1994) 

54 P. Markos: J. Phys.: Condens. Matt. 7, 8361 (1995) 



G 



